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Droplet patterns condensing on solid substrates (breath figures) tend to evolve into a self-similar regime, char¬ 
acterized by a bimodal droplet size distribution. The distributions comprise a bell-shaped peak of monodisperse 
large droplets, and a broad range of smaller droplets. The size distribution of the latter follows a scaling law 
characterized by a non-trivial polydispersity exponent. We present here a numerical model for three-dimensional 
droplets on a one-dimensional substrate (fiber) that accounts for droplet nucleation, growth and merging. The 
polydispersity exponent retrieved using this model is not universal. Rather it depends on the microscopic details 
of droplet nucleation and merging. In addition, its values consistently differ from the theoretical prediction by 
Blackman (Phys. Rev. Lett., 2000). Possible causes of this discrepancy are pointed out. 

PACS numbers: 68.43.Jk, 47.55.D-, 89.75.Da, 05.65.+b 


I. INTRODUCTION 

When a flux of supersaturated vapor gets into contact with a 
solid substrate, a condensation process can originate, leading 
to the formation of droplets patterns on the substrate (“breath 
flgures” ID). The interest for breath flgures is both theoreti¬ 
cal and practical. From the theoretical point of view, they can 
be used as a test ground for scaling concepts in a well-posed 
non-equilibrium setting. From the practical point of view, they 
appear in many natural phenomena, e.g. dew deposition on a 
spider net or on a leaf. They can also be exploited for techno¬ 
logical applications: water collection from dew harvesting p2L 
m, biological sterilization Q, manufacturing of surface struc¬ 
tures and patterns for nano-technologies fabrication of 

efficient heat-exchangers and cooling devices ifTOlfT^ . For 
such applications, understanding and controlling the droplets 
formation, growth and coalescence is crucial. 

The formation of breath flgures undergoes several phases 
ca. First, the droplets nucleate on the substrate; then they 
grow and coalesce, creating a roughly monodisperse distribu¬ 
tion. Eventually the space released by merging is sufficient 
for the nucleation of new droplets. As the evolution contin¬ 
ues, self-similar droplet patterns appear EHni. 

In this phase, bimodal droplet size distributions emerge 
in experiments |[T4l fT8ti2T1l . as well as in simulations 1211 - 
1^ . The size distributions feature a monodisperse bell-shaped 
peak for the largest droplets, and a power-law distribution for 
the smaller droplets, characterized by a non-trivial polydis¬ 
persity exponent. Scaling descriptions for the droplet num¬ 
ber density have been largely adopted in the classical theory 
for breath flgures CH [15] |22l [24H271 . Such a theory relates 
the polydispersity exponent to the exponents for the time de¬ 
cay of the droplet number and the porosity, i.e. the fraction 
of the non-wetted area over the total area of the substrate 
csiini. The scaling descriptions of the droplet size distri¬ 
bution are also solutions (281 [29l of the Smoluchowski coagu¬ 
lation equation (301, an integro-differential equation describ¬ 
ing the evolution of the droplet size distribution as a conse¬ 
quence of merging processes. However, for droplets grow¬ 
ing on two-dimensional substrates the polydispersity expo¬ 
nent derived from the solution of the Smoluchowski coagu¬ 


lation equation ED was found to be noticeably smaller than 
the value found in numerical simulations EDEl and larger 
than the one observed in experiments (211 . 

In the present paper, we examine the case of a one¬ 
dimensional substrate with three-dimensional droplets. We 
use a numerical approach in order to test the existing scal¬ 
ing theory and, in particular, the common assumption that 
the polydispersity exponent takes a universal value (Hi [T5l 
|22l[24H27l|3ll. We introduce a numerical model based on a 
nucleation rule governed by surface tension-driven instabili¬ 
ties. The dynamics account for the presence of a precursor 
film between droplets, droplet growth due both to direct mass 
deposition from the surrounding vapor and surface diffusion, 
and non-trivial droplet interactions due to deviations from the 
spherical shape. By means of simulations, we systematically 
explore the dependence of the droplet patterns on the depo¬ 
sition rate, the rules of droplet interaction, the radius of the 
flber and the nucleation radius. Surprisingly, we And a de¬ 
pendence of the polydispersity exponent of the droplet size 
distribution on the microscopic details of the model, leading 
to the conclusion that the polydispersity exponent is not uni¬ 
versal, as assumed in the classical scaling models. Moreover, 
we observe a sizable mismatch between the predicted ED and 
the observed values of polydispersity exponent. We point out 
possible sources of the discrepancies. 

The paper is organized as follows. In Sect. |n| we intro 


duce our model and its numerical implementation. In Sect. Ill 


we revisit the scaling description of breath flgures with spe¬ 
cial emphasis on the predictions for droplet growth on flbers. 
These predictions are then, in Sect. |IV| carefully tested by 
comparison to comprehensive numerical data. Finally, we 
conclude in Sect. [Vl 


II. MODEL AND NUMERICAL METHOD 

The precise dynamics of droplet nucleation and growth on 
a solid substrate is still a debated matter and depends on the 
specifle system. Possible mechanisms include Eda hetero¬ 
geneous nucleation on impurities and defects of the substrate, 
homogeneous nucleation, hydrodynamic instabilities, surface 
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FIG. 1. Sketch of the model adopted for the numerical simulation of 
the growth of droplets % and z + 1 with radii Ri and Ri+i, respec¬ 
tively. The droplets reside on a fiber (solid gray line). They grow due 
to a constant water flux per unit length, $ [/xrn^s”^], due to the su¬ 
persaturated vapor impinging on the fiber. There is a precursor film 
of radial thickness Ci in the region between adjacent droplets, which 
has a length li. When the precursor film collects sufficient water vol¬ 
ume, new droplets are nucleated. When two droplets approach each 
other within their interaction ranges sRi and sRi-^i, they merge. 

diffusion on the substrate, direct condensation of the vapor 
molecules on the droplets surface, coalescence with neigh¬ 
boring droplets, heat transfer localized at the triple line and 
preventing nucleation. The radius of an isolated droplet on a 
substrate increases following a power law in time, whose ex¬ 
ponent changes according to the specific growth mechanism 
(1311331 [34l. However, in densely populated droplets systems, 
the growth is always dominated by coalescence. Therefore, 
the growth exponent of the average radius depends only on 
the dimensionalities of the droplets and the substrate da. 

A. Modeling breath figures on a fiber 

We model breath figures on a fiber as a one-dimensional 
chain of spherical droplets, in contact with a supersaturated 
vapor. The droplet has a radius Ri and an interaction range 
sRi, where 6 is constant for all the droplets. Such an interac¬ 
tion range has been introduced to keep into account the devi¬ 
ation from the spherical shape observed in experiments Ea 
(see Fig. [^t. We account for a cylindrical water prewetting 
film surrounding the fiber, in the spaces between the droplets. 
Water is deposited on the fiber by condensation from the sur¬ 
rounding vapor. We do not consider removal of water from 
the fiber by evaporation and gravity. 

We consider a case where the supersaturation of the va¬ 
por is small and the diffusive transport to the thread is the 
growth limiting process for the droplets. In this case we ex¬ 
pect cylindrically symmetric vapor concentration profiles in 
the far field. We assume that the system is in thermal steady- 
state conditions, i.e. the water flux deposited on the fiber T> - 
water volume per unit length per unit time - is small enough to 
keep the temperature of the droplets and the fiber uniform and 
invariant in time. We also assume hydrodynamic steady-state 
conditions, taking T> to be constant in time and independent of 
the position, x, on the fiber. These approximations have been 
largely adopted in the study of breath figures ifT^ and can be 
considered reasonable for laminar vapor flows. The modelling 


of more complicated vapor flows may require the inclusion of 
correction terms, but this is beyond the scope of the present 
paper. 

B. Evolution of the droplets 

We adopt an event-driven approach and periodic boundary 
conditions for the flber. In order to simplify the problem, we 
decouple the treatment of the droplets growth due to the im¬ 
pinging flux from the nucleation process. In particular, for 
each time step, we consider the growth in a time-continuous 
fashion and the nucleation in a time-discrete fashion. The 
algorithm to advance from time tj to the next instant fj+i 
proceeds as follows. At first, we only consider the droplets 
growth due to the water flux impinging on the droplets them¬ 
selves and we disregard the nucleation of new droplets. The 
constant flux density T> impinging on a length 2Ri covered on 
the flber by the particle i results in a growth law 

s 

We integrate Eq. Q in time, from ttot At, finding 

I ^ At 

+ = ( 2 ) 

We calculate the time intervals Ati for binary merging events 
between adjacent droplets i and i + 1, by solving the system 

— Xi = {1 s) [Ri-^i{t -b Ati) + Ri{t + ^ti)] ( 3 ) 

where Xi and are the positions of the respective cen¬ 
ters (note that they only move when the droplets merge). We 
take the minimum of these time intervals, Atmerge = niin 
as a first estimate of the time step to advance the system in 
time. This time step corresponds to the first merging event 
that would take place if there was no nucleation and no water 
deposited between the droplets. We calculate the volume of 
water V"gap,i deposited on the gaps between adjacent droplets i 
and i + 1 during the time Atmerge, by integrating 

= ^[Xi+i -Xi- Ri+i - Ri] ( 4 ) 

between t and t + Atmerge- 

C. Nucleation of droplets 

In each time step we determine if one or more nucleation 
events could happen during the time interval Atmerge- To 
this aim, we consider the existence of a cylindrical precur¬ 
sor film, surrounding the parts of the flber where no droplets 
are present. The precursor film grows in time by deposition 
of mass and may eventually develop into a surface tension- 
driven instability. Stability analysis ll^ for a flber of radius 
Rf reveals that perturbations with a wavelength A > 27tR f 
are unstable. The most unstable perturbation, i.e. the fastest 
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FIG. 2. (Color online.) Time evolution of the space occupied by the droplets for e = 2%, T> = lOOO/xm^/s, Nq = 1.5 x 10^, and 
Rf = Rmin = 1 A^ni. The different colors represent the different ages of the droplets at the considered time instant and they vary from dark 
(blue online), for the youngest and still very small droplets to lit-gray (yellow, then red online), for the oldest hence largest droplets. The 
central panel is a magnification of the area indicated by the black frame in the left panel, and the right panel shows a magnification of the frame 
in the central panel. They are both displayed on a logarithmic scale with the origin corresponding to the instant when the free area started to 
be populated {Le. the bottom side of the black frame in the previous panel). 


growing perturbation, has a wavelength A* = 2^/27tR f ll36l . 
Therefore, A* can be regarded as a characteristic wavelength 
of the system. Necklace-shaped chains of droplets have been 
observed in experiments on fibers, where equispaced droplets 
appeared at a distance A* from each other ||36l. In our model, 
we consider that a nucleation event can take place between 
two droplets i and i + 1, during the time Atmerge only if the 
gap length li between such droplets, is larger than the charac¬ 
teristic wavelength of the first unstable perturbation A*. The 
gap length is defined as(t) = x — 

and, in order to have a nucleation event, it has to be 

li> X* = 2V2TrRf (5) 

From the length of the gap we determine the number of 
possible equispaced nucleation sites in the gap, Ni = [Z^/A* J. 
Additionally, we impose a minimum droplet size Rmin^ for 
the nucleation event to take place. We approximate the film 
between the droplets i and i -j- 1 as a cylinder with volume 

Vfilm.i = TT [{Rf + Cif - Rj] k , (6) 

where is the thickness of the film. We check if the amount 
of liquid V"gap,i deposited in the gap li during the time Atmerge 
is large enough to satisfy the following condition: 

+ fgap,i ^ 5 

where ^ is the volume of the precursor film in the gap 
at the beginning of the time step Atmerge. and V^ap,i is derived 
by integrating Eq. 0 over Atmerge- In nucleation processes, 
the critical droplet radius required to have a nucleation site 
can be calculated as ISTll Rnud = 2cr/ [plRvT^^{pv/P oo)], 
in which a is the surface tension between the liquid and the 
vapor, pl is the density of the liquid, T is the temperature ex¬ 
pressed in Kelvin, Ry is the universal gas constant for water 
vapor, py /poo is the supersaturation rate, withpy the pressure 


of the vapor and poo the pressure of the saturated vapor. In our 
model we adopt the minimum droplet size as i?min = ^nuci- 
For typical environmental temperatures (20-40°C) and su¬ 
persaturation rates py /poo > 1.001 we have Rmxd < 1 pni. If 
Eq. Q and Eq. (|7]) are both satisfied, a new nucleation event 
takes place. If Eq. ^ is satisfied but Eq. (|^ is not, we con¬ 
sider the volume of water lgap,i deposited on the gap as 
increasing the thickness of the cylindrical precursor film, and 
we store the information until the next time step. If neither 
Eq. (0 nor Eq. (|7]) are satisfied, we consider the amount of 
water deposited on the gap Vgap,i as collected by the droplets 
adjacent to the bridge itself. In particular, this water volume 
will be collected by the largest droplet, due to the Laplace 
pressure pl EH. Such a pressure depends on the local cur¬ 
vature and, in the case of a spherical droplet of radius R, it is 
given by Pl = 2afR. 


D. Merging of droplets 

We advance the system in time, by evolving it accordingly 
to the estimated Atmerge- we grow the droplets, we insert the 
new nucleated ones and we merge the droplets when they ap¬ 
proach each other within their interaction ranges. In order 
to do so, we replace the two merging droplets with a new 
one, with a mass equal to the sum of the masses of the merg¬ 
ing ones. Its center is located in the center of mass of the 
two droplets. In the simulation we explicitly prevent triple 
merging events between neighboring droplets. Though rare, 
such unphysical events could occur in the described numeri¬ 
cal scheme, due to the collection of liquid from the adjacent 
water bridges. When triple merging is at hand at the end of 
the estimated time step Atmerge. we halve the estimated time 
step, and we evolve the system accordingly. 
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FIG. 3. Droplet number density for £ = 2%, Rf = Rmin = 1 Mm, 
4> = lODO/im^/s, Lf = I m. No = 1.5 x 10^. The solid black line 
represents n(s, t). The dashed black line is a fitting line for the poly- 
disperse scaling range of the droplets and it has a slope — r on the 
log-log scale. The gray solid line represents the cutoff for the large 
droplets /(s/E(t)). The gray dashed line represents the cutoff for 
the small droplets: g{s/so). The two vertical dash-dotted lines de¬ 
limit the scaling range. For small droplets (s < A*) the distribution 
is dominated by g{s/so). For large droplets (s/E(t) > 1%), the dis¬ 
tribution is dominated by /(s/E(t)). In this range, we can identify 
a monodisperse bump, corresponding to the oldest droplets, and a 
gap, corresponding to the times where the spaces released by merg¬ 
ing droplets were not large enough to allow the nucleation of new 
droplets. The present graphs are derived by averaging the droplet size 
distributions over 10 time instants and over five simulations with dif¬ 
ferent but equivalent initial conditions. The considered instants are 
in the self-similar regime (t > 1000 s) and they are chosen in such a 
way to have 100 points per time decade. 
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FIG. 4. Sketch of the simplified model adopted for deriving the time 
growth law of the size of the largest droplet E(t). The exponent of 
such a time growth is calculated based on the assumptions that the 
droplets are monodisperse and in touch with each other. The consid¬ 
ered mechanisms governing the growth are the droplet coalescence 
two by two and the deposition of water from a uniform external fiux 
per unit length 4>. 


scale cutoff accounts for the minimum droplet size i^min and 
the characteristic length of the surface tension-driven instabil¬ 
ity in the precursor film A*. The large-scale cutoff comprises 
a bump and a gap (see Fig.[^. The bump represents the oldest 
droplets in the system (first droplet population). The gap orig¬ 
inates from the times where the openings generated between 
the first droplets were still too small to admit a second wave 
of nucleation. 


III. SCALING THEORY FOR DROPLETS ON A FIBER 

In Fig. we show the time evolution of the droplets and 
their coverage of the fiber surface. The horizontal colored 
segments represent the areas covered by the droplets on the 
fiber at a certain time t. Different colors reflect the differ¬ 
ent ages of the droplets at the speciflc time t. In black (blue 
online) we represent newly created small droplets, and the 
color fades (turning to yellow, then red online) as the droplets 
become older. When two droplets merge, they release two 

regions of length ^ Ri R 2 — {Ri + R^)^ ^. For suffi¬ 
ciently big droplets the released space is large enough to host 
new droplets. The opening gaps become larger with the in¬ 
crease of the largest droplets. Eventually, new droplets nu¬ 
cleate and grow in the gaps, forming a hierarchical struc¬ 
ture in the droplet size distribution. After a sufficiently long 
time, self-similar droplet patterns emerge: pictures taken at 
different time instants look alike, apart from length rescaling 
(c/Fig.|^. This is reflected in a power-law size distribution 
of the middle-sized droplets |[T5l[^[24H27l . characterized by 
an exponent r (see Fig. [^, namely the polydispersity expo¬ 
nent. On the other hand, the largest (oldest) and the smallest 
(newest) droplets have a different behavior. Their physics is 
captured in cutoff functions which describe the termination of 
the power law for large and small droplets. 1^ . The small- 


A. Size of the largest droplets 

We consider the following simplifled setting, in order to 
give an estimate of the growth law in time of the largest 
droplet size E(t). At time we take a chain of monodisperse 
droplets of radius Rn and size E(f^) = Rf^, in contact with 
each other (see Fig.|^. At time t*, the droplets have merged 
two by two, therefore E(t*) = 2E(tn)- Since the merging 
process is almost instantaneous, t* After merging, the 

droplets start to grow as an effect of the deposition of a uni¬ 
form water flux per unit length, on the flber. Due to mass 
conservation, at times tn and tn+i, it must be 2R^t = ^ttR^. 
Hence, 

/ 3 ^\ 

E=f^j with ^ = 3/2. (8) 

The so-calculated exponent 2 : is in agreement with both ex¬ 
perimental flndings (SEl and theoretical predictions ifBl . The 

prefactor (|f has to be regarded as a lower bound. Since 
the merging is assumed to happen instantaneously, the largest 
droplet in the system at time t* will most likely be a 

droplet that has just originated from a merging. Hence, the 

upper bound for the prefactor will be 2 (|f 
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B. The droplet number density n{s,t) 

The droplet number density n{s,t) is defined as the number 
of droplets of size s at time t per size interval ds of droplets 
and unit length of the substrate, where s = is the size 
of a droplet with radius r. In our case, n{s,t) has the units 
of Consequently, the Buckingham-Pi theorem |[39l 

ensures that the droplet number density can be expressed as 
n{s,t) = s~^ s/sq), where f{x) is a dimension¬ 

less scaling function, sq characterizes the cutoff for the small 
droplets, 'E{t) is the maximum droplet size at time f, and 
0 = {D d)/D 1221 is an exponent depending on the di¬ 
mensionality of the droplets D and of the substrate d. Here, 
we consider the case of three-dimensional droplets, D = 3, 
growing on a one-dimensional fiber, d = 1. Therefore, from 
merely dimensional considerations one infers that 0 = 4/3. 

The classical scaling theory for breath figures (151 ESI es¬ 
tablished that, in the late regime, the droplet size distribution 
becomes self-similar: there is an increasing scaling range be¬ 
tween So and S(t), characterized by a polydispersity exponent 
r (Fig.[^. This suggests that the function f{s/T^{t), s/sq) can 
be factorized into a power law and two cutoff 

functions, f {s/T^{t)) and ^ (s/sq), accounting for the large 
and the small scales, respectively. For their asymptotics we 
request that f{x) = /o =const for x ^ 1, and that g{x) = 1 
for X ^ 1. 

Thus, the droplet number density is expressed as 


n(s, t) 


c-o 


m 


e-r 




(9) 


In the framework of the renormalization group theory 
(401 l4T]| . one would expect r to be a universal constant, de¬ 
pending only on the dimensionality of the system and not on 
microscopic details The value of r is 

related to the decaying exponents characterizing the time evo¬ 
lution of the porosity and the number of droplets, by consis¬ 
tency reasons (see Sections III C| and IIID| ). This poses some 
limitations to the range of the physically acceptable values of 
r. Within such range, a theoretical derivation of r has been 
proposed by Blackman and Brochard (STI . In Section [HT eI 
we will revisit it, then we will compare the expected value of 
r to our numerical results (Section |IV]). 


a constant /o otherwise. Similarly, the role of g{x) is to pro¬ 
vide a cutoff for small droplet sizes. Therefore, it takes the 
value g{x) = 1 for all but the smallest values of x. Conse¬ 
quently, 

pXN 

N{t) = fo / x-^dx 

d Sl\ 




fo ^1 — 0 


1 — r 


1-^ _ f sn a 


1 — 7 




( 11 ) 


where x at — 1 and sn — sq are constant. Substituting Eq. ^ 
into Eq. ^TT\ we find 


N{t) - 


fo 


- 

Xn 






im 


V ^ 


( 12 ) 

and fort^oc the number of droplets hence decays in time 
as 


N{t) - r 


with an exponent 


^/^U(('-l), ifT<l, 
1 z{9 — t ), if t > 1 . 


(13a) 


(13b) 


We note that the values of the exponent z' do not depend 
on the specific choice of sn, provided that there is a suffi¬ 
cient scale separation between x^ and The case of 

r < 1 corresponds to a monodisperse droplet population. The 
corresponding trivial exponent z' = 1/2 can also be found 
from the consideration that an ideally monodisperse popula¬ 
tion of droplets of size S will cover the whole length of the 
fiber Lf, in the limit t ^ oo. Therefore, in such a case, one 
would have 1 ^ ^ such that N ^ . In 

particular, this exponent describes the decay of the number of 
large droplets populating the bump of the roughly monodis¬ 
perse large droplets in the large-scale cutoff function (Fi g. [3]) . 
For a poly disperse droplet population r > 1, and Eq. ( |1%[ ) 
sets an upper limit for the range of the physically acceptable 
values of r. The size of the smallest droplets is fixed and 
the larger droplets grow. Hence, the total number of droplets 
must decay and z' > fi (see Eqs. (13)). Since > 0, it must 
be 1 < r < 6>. 


C. The number of droplets per unit length, N (t) 

The total number of droplets N{t) at time t per unit length 
of fiber can be written as 


N{t) = r 

Jo 


n{s, t) ds 


( 10 ) 


To evaluate the integral we substitute Eq. ([^ into Eq. ( p^ , 
and we define a new variable x = We replace the 

function inside the integral with a proper combination of step 
functions, using the fact that f{x) should contribute to the 
shape of the droplet distribution only for large droplets, being 


D. The porosity p{t) 

The porosity p is defined as p = 1 — AdjAtou where Ad 
is the wetted area covered by the droplets and is the total 
area of the substrate. For one-dimensional fibers these quan¬ 
tities correspond to the wetted length and the total length of 
the fiber, respectively. It is convenient to define an effective 
porosity p* that keeps into account the interaction ranges of 
the droplets 


N 


— 1 — Ai{l -h s)/At^ 


(14) 


i=l 
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Here, Ai {1 e) is the effective area occupied by the 
droplet, i.e. for a one-dimensional substrate, Ai = 2Ri, 
with Ri the radius and (eRi) the interaction range of the 
droplet. In the general case Ai ^ such that the effective 
porosity is 

poo 

p*(t) ^ 1 — / n(s, t) (1 + 5 ) ds , (15) 

Jo 

where C is a constant depending on the geometry of the sys¬ 
tem. Following a procedure similar to the one adopted to cal¬ 
culate the number of droplets per unit length N{t), and intro¬ 
ducing the exponent 0 = {d-\- D)/D,wq find 

p*it) = l--^C{l + e) f ” dx, 

^ Jsp/^it) 

where again Xp 1 and Sp ^ sq, even though they may 
slightly differ from the integration limits used to calculate the 
total number of droplets, xn and s^. After few algebraic 
steps, we derive 


p*{t) = 1 


fo 

e-T 


X 


e-T 

p 


s 


e-T 

p 




C(l+e). (17) 


In the late regime, t ^ 00 , we expect p*(t) ^ 0, because 
the number of gaps decays together with the droplet number, 
Eq. ( |13b| ), and the length of the gaps is bounded in our system. 
Hence, 1 - /o C (1 + 5 ) - r) = 0 and 


p*{t) c. 


/o 4-" 




ZTY c {I+ e), 


z{e-T) 


( 18 ) 


Therefore, in the late-time scaling regime, the effective poros¬ 
ity decays in time as 


p%t)r^t-^ withk = z{0-r). (19) 


Note that this exponent takes the same value as the exponent z' 
determined for the N{t), when r > 1 {cf. Eqs. @). This in¬ 
dicates that the porosity can be viewed as the number of gaps 
(coinciding with the number of droplets for a one-dimensional 
system) multiplied by a characteristic gap size, which takes a 
constant value in the long-time limit. 

Interestingly, this result also holds when considering only 
droplets larger than a fixed finite size Sc > sq, i.e. smaller 
droplets are considered part of the larger gaps. One can check 
by straightforward calculations that only the prefactors of the 
asymptotic power laws are affected by this change of the lower 
cutoff. 


E. Expected value of r 

The value of r has been theoretically determined by Black¬ 
man and Brochard ED from an analysis of the scaling of 
Smoluchovski’s equation |[3Ql for droplet coagulation, which 
is then evaluated in a renormalization group framework under 


the assumption that the particle number and the porosity must 
obey the same time dependence. 

The derivation of Blackman and Brochard ED relies on the 
assumption that the collision rate ncoii(5i, S 2 ; t), i.e. the num¬ 
ber density of droplets of size si and S 2 that collide at time 
t per unit length per unit time, factorizes into the product of 
single particle distribution functions and a geometrical factor: 

ncoii(si,S2;i) ~ (20) 

where ncoii(<si, S 2 ; t) is the number density of collisions per 
unit length of substrate between two droplets of size si and 
52, from time 0 to time f; n{siA) is the number density of 
droplets of size si per unit length of substrate at time t and it 
represents the probability density of having a droplet of size 
5i at time t (the normalization constant Lf can be seen as 
equivalent to the number of knots in a Boltzmann lattice); 
n(52, t) / A^(t) is the probability density of having a neighbor¬ 
ing droplet of size 52 , uncorrelated to the size si. The term 
represents the speed at which the edges of 
the two adjacent droplets approach each other, by effect of 
the mass deposition through In order to demonstrate this, 
in the case of D = 3 and d = 1, one can write such speed 
as —/i, where the gap between the two adjacent droplets is 
h = X 2 — xi — {R 2 + ^i)(l + ^)- By substituting Eq. ^ 
and taking the derivative with respect to time, one then finds 
Ri - 

The total collision rate per unit length of substrate at time t, 
Nco\\{t) = f^ 'dco\\{si, SjA) dsi dsj is expected to scale 
as ISTl A^coii(^) ^ ^(^)- equating such a scaling as¬ 
sumption and Eq. ( [20| ), upon substitution of Eq. ( p3] ), a re¬ 
lationship is found, among r, d and D. A similar procedure 
is repeated, keeping into account Acoii(f), the change in the 
substrate area covered by the droplets per unit length, due to 
collision events, instead of Acoii(^)- The combination of the 
relationships among r, d and D obtained through such a scal¬ 
ing analysis yields ED that Ttheor = 7/6 for the growth of 
three-dimensional droplets on a fiber. We refer the reader to 
ED for the details of the derivation. 


IV. RESULTS 

A. Setup of simulations and growth of the largest droplets 

All simulations start from an initial condition where Nq = 
1.5 X 10^ droplets form an equispaced necklace at the dis¬ 
tance A* from each other; the length of the fiber T j is adapted 
accordingly. The initial conditions differ by a slight poly- 
dispersity of the droplet radii, which are chosen as ro,i = 
Rmin{f + 0.01 /rand,i) where I < i < Nq, and /rand,i are ran¬ 
dom numbers in the interval 0 < /rand ,2 < 1 • We select the 
characteristic size of the distribution, E(f), to be the maxi¬ 
mum droplet size at time t. As predicted by Eq. it scales 
as E ~ (If ^ ~ g. The asymptotic value 

of 2 in the reduced plot indicates that most likely the largest 
droplet in the system has recently undergone a collision. 
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FIG. 5. Time evolution of the maximum droplet size (same parame¬ 
ters as in Fig.[^. The black solid line represents the values inferred 
from our simulations. The black dashed line has a slope of 3/2 on 
the double logarithmic axes, proving that E ~ and it represents 
the function y — (3T>t/27r)^. The gray solid line is a reduced plot, 
obtained by dividing E by the RHS (right hand side) of Eq. j8}. 


FIG. 6. Time evolution of the effective porosity (same parameters 
as in Fig. [^. The lower lines (black) represent the total effective 
porosity, as defined in Eq. (Td] ), for five simulations with different but 
equivalent initial conditions; the upper lines (gray, solid) represent 
the partial effective porosity p* — 1 — A*/A^ou where A* is the sum 
of the areas (lengths) covered by the droplets of size s > A*^. The 
gray dashed line represent the line used to fit the porosity p* and it 
has a slope —k = —0.332 zb 0.003, on the log-log axes. Erom this 
we estimate r = 1.112 zb 0.002. 


B. Scaling of the droplet number and the porosity 


We check the consistency of the expon ents for p* (t) and 
N(t) by inspecting Eq. ( p^ and Eq. (13b), respectively. 

In order to improve the statistics, we run five simulations for 
each case, with different but equivalent initial conditions. To 
this aim, we take different seeds in the random number gen¬ 
erator used to create the initial distribution of the radii. The 
curves p*(t) and N(t) for different random numbers seeds 
all lie on top of each other, and they have the same scatter¬ 
ing of the data (see Eigs. and|7]). The five runs produced 
data at time instants which were not perfectly in-phase, due to 
the event-driven nature of the model. Therefore, we use the 
curves p*(t) and N(t), produced by overlapping the individ¬ 
ual curves of the five runs, to extract the exponents k and z' 
respectively. 

In Eig.j^we show the time evolution of the modified poros¬ 
ity p*(t), for five equivalent initial conditions. The following 
parameters were used: 5 = 2%, Rf = Rmm = 1 ^ = 

1000/im^/s. The lower lines (black) represent the total ef¬ 
fective porosity, as defined in Eq. The upper lines (gray, 
solid) represent the partial effective porosity p* = 1—^4* /^Ttot, 
where ^4* is the sum of the areas (lengths) covered by the 
droplets of middle and large size. Specifically, we keep into 
account only the droplets of size Sc > A*^ and we use these 
curves to find the exponent k. As pointed out at the end of 
Sect. HID this choice does not affect the resulting exponent, 
but it reduces the fluctuations related to the microscopic de¬ 
tails of the nucleation of new droplets. The gray dashed line 
shows the best fit of the slope —k calculated from fitting the 
modified porosities p* of the five simulations. Eor the consid¬ 
ered parameters, we obtain k = 0.332 ± 0.003. 

In Eig. [T] we show the time evolution of the number of 
droplets, for the same data of Eig. The dark gray thin 
lines (blue online) represent the total number of droplets N{t). 
The black lines represent the number of the large droplets: 
s/E(f) > 0.75^; by fitting them we estimate a decaying ex¬ 


ponent of 0.503 ± 0.026 ^ 1/2, in agreement with the theo¬ 
retical prediction for a monodisperse droplet populations (see 
Eq. ( [13b| ), for r < 1). The lit-gray thin lines (green online) de¬ 
pict the number of the small droplets, with s/'E{t) < XNo{t) 
and XNo{t) = A*^/E(t). They suffer from large oscillations, 
reflecting the repopulation of large areas that are episodically 
released by the merging of large droplets. The gray thick lines 
(red online) show the number of droplets of middle and large 
size, with s > A*^, as defined above. We use these lines to 
And the exponent characterizing the decay of the number of 
droplets in time z' = 0.329 ± 0.002. Erom Eq. (|13b| and 
Eq. we expect k = z', which is verified within the esti¬ 
mated error; the relative error on the average value is ^ 0.4%. 
Erom Eq. we derive the polydispersity exponent of the 
droplet size distribution r = 1.112 zb 0 .002, which is within 
the physically acceptable range 1 < r < 0. 


C. Droplet number density 

We calculate n{sR) from our numerical data by dividing 
the range of the sizes s into Abins bins of width Asj, cen¬ 
tered around Sj, where j = l,..,Abins. At each time in¬ 
stant, we then calculate the droplet number density n{sjR) 
by counting how many droplets lie in the respective bins, 
{sj — Asj/2) < s* < {sj -f Asj/2), and dividing the re¬ 
spective numbers by (LfAsj), where Lf is the total length 
of the fiber. We then take the average of the resulting droplet 
size distributions over 10 time instants and over five simula¬ 
tions with different but equivalent initial conditions, each one 
with initial droplet number Nq = 1.5 x 10 ^. The considered 
time instants are in the self-similar regime (t > 1000 s) and 
they have been chosen in such a way to have 100 points per 
time decade. In particular they belong to the time range 1000 
s<t< 2000 s. 
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FIG. 7. (Color online) Time evolution of the number of droplets 
per unit length (same parameters as in Fig. [^. The dark gray thin 
lines (blue online, uppermost right, denoted as ’’TOTAL” in the leg¬ 
end) represent the total number of droplets iV(t), for 5 simulations 
with different but equivalent initial conditions. The lit-gray thin lines 
(green online) represent the number of the small droplets, with s < 
and its large oscillations. The mid-gray thick lines (red online) 
represent the number of large and middle size droplets, with s > ; 

upon fitting, we derive a decaying exponent z' = 0.329 zb 0.002. 
The black lines represent the number of the large droplets, with 
s/E(t) > 0.75^; their decaying exponent is 0.503zb0.026 ~ 1/2, in 
agreement with the theoretical prediction for a monodisperse droplet 
population. 



FIG. 8. Droplet number density n{s,t) (same parameters as in 
Fig. 1^. The black line represents n{s,t)s^. The gray solid line 
represents the same quantity rescaled: n(s, t)s^(s/E)'^“^, with 
T = 1.112, from the porosity fit: such an estimate seems to be 
consistent, since the plotted curve is horizontal, in the self-similar 
range (middle size droplets). The gray dashed line represents the 
same rescaled quantity with Ttheor = 7/6 ISTl : such a value is not 
consistent, since the plotted curve is not horizontal, in the range of 
self-similar sizes. The present graphs are derived by averaging the 
droplet size distributions over 10 time instants and over five simula¬ 
tions with different but equivalent initial conditions. The considered 
instants are in the self-similar regime (t > 1000s) and they are cho¬ 
sen in such a way to have 100 points per time decade. 


The resulting distribution n(s, f) is shown in Fig. (black 
solid line). The gray solid line and the gray dashed line repre¬ 
sent our numerical estimate of the cutoffs for the large droplets 
/(s/E(t)) and for the small droplets g{s/so), respectively. 
They are obtained from Eq. (|^, with the specific choice for the 
cutoff asymptotes /(s/E(f)) = /o =const for s/E(t) < 1%, 
and g{s/so) = 1 for s > A*. 

In Fig. we show the droplet number density rescaled in 
such a way to make it dimensionless. The black solid line 
represents the quantity n{s,t)s^. With the gray solid line we 
introduce a further rescaling n{s^t)s^{s/'E{t)y~^, in order 
to verify the estimate of r = 1.112 from the fit of the effec¬ 
tive porosity. The gray dashed line depicts the same quantity, 
rescaled with the theoretical prediction Ttheor = 7/6 lISTIl . The 
gray solid line is horizontal, in the range of sizes correspond¬ 
ing to the self-similar regime, while the gray dashed line has a 
slight positive slope. Consistently, our data support a polydis- 
persity exponent considerably smaller than expected, but still 
matching the values obtained from the decay of the porosity 
Fig. and the droplet number Fig. [ 7 ] The good agreement 
fully supports the classical scaling results relating the power- 
law dependencies of these quantities 


D. Parameter dependence of r 

We repeated the procedure to calculate r and to verify its 
accuracy, varying several parameters, in order to clarify its 
nature in terms of universality. In particular, we considered 


its dependence on the interaction range 5 between droplets, 
the nucleation radius i^nuch the fiber radius Rf, the imping¬ 
ing water flux per unit length T> and the monodispersity of the 
nucleated droplets chain. The variation of the latter, has been 
realized by randomly redistributing different percentages of 
the water volumes V^ap,i accumulated on the gap, among 
the new nucleated droplets on the same gap. We found that 
changes in both the impinging flux T> and the monodispersity 
of the nucleated droplets do not affect the exponent r, in line 
with expectations (see Sect. However, the exponent r ex¬ 
hibits a clear dependence on the interaction range, as it grows 
with 5 (see Fig.j^a), solid line). The ratio between the nucle¬ 
ation radius Rmci and the characteristic length of the perturba¬ 
tion A* also has an influence on r. In particular, r decreases 
with increasing i7nuci/A* (see Fig. ib), solid line). 

Remarkably, all measured values of r are substantially dif¬ 
ferent (4%-8%) from the theoretical prediction ED demand- 
ing Ttheor = 7/6 (see Fig. dashed lines). The difference is 
signiflcant, since it exceeds 10 times the estimated error of the 
exponent. 

The dependence of the r exponent on the interaction range 
5 can be explained as follows. Having larger interaction 
ranges 6 implies having a larger area released upon merg¬ 
ing. One can see this by considering two merging droplets 
of radius Ri and 7^2, with interaction range sRi and eR 2 
respectively. Just before the merging, they cover an area 
A = 2(Ri + 7^2)(1 + 5 ). After the merging the covered area 
will be A' = 2(1 + 5 ) 77 ', where 77' is the radius of the new 
generated droplet. Therefore, the released area, i.e. the gen¬ 
erated gap, will be Ag = 2(1 + 5)(77i + 772 — 77'). Larger 
generated gaps imply an enhancement in the nucleation pro- 
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FIG. 9. The exponent r as a function of (a) the interaction range, £, and (b) the ratio among the radius of the nucleated droplets Rnud and 
the characteristic length of the perturbation originating the droplets chain A*. The solid lines represent the values of r inferred from the fit 
of the effective porosity. The dashed line is the theoretical value predicted by the classical theory for breath figures |3ll Ttheor = 7/6. The 
parameters in (a) are Rf = i^min = 1 A^m, T> = 1000 /im^/s, and iVo = 1.5 x 10^, while for (b) we used £ = 0%, $ = 1000 /xm^/s, and 
No = 1.5 X 10^ 


cess, because more space and more water volume are available 
for the new droplets. On the other hand larger gaps imply a 
delay in further merging processes. Such effects are particu¬ 
larly relevant when the merging droplets are large, i.e. when 
{Ri R 2 — R') is large. For increasing values of 5 , one can 
then expect an increase in the number of small droplets and a 
decrease in the number of large droplets. Since — r represents 
the slope of the size distribution of the droplets n{sR) (see 
Fig. [^, this will result into higher values of r (see Fig. |^a), 
solid line). 

The dependence of r on the ratio i^nuci/^* can be explained 
through the following qualitative argument. We consider the 
case where the nucleation radius Rnuci changes, but the char¬ 
acteristic length of the perturbation A* remains constant, as 
well as all the other parameters. Given a certain newly gen¬ 
erated gap of size the number of possible nucleation sites 
Ni inside such a gap does not change with Rnuch but only 
with A*, since Ni = The time interval At during 

which Ni droplets can nucleate is the time required to reduce 
the gap size from to A^^A*. After that time, only Ni — 1 
droplets will have the space to nucleate inside the gap. How¬ 
ever, the water volume required to have a nucleating chain of 
Ni droplet is larger for larger i^nuci and criterion 0 may not 
be matched during the time interval At. Therefore, for larger 
Rnud, a smaller number of new droplets will appear, before the 


gap closes. Overall, a lower number of small droplets will be 
present in the system, hence a larger number of big droplets, 
due to mass conservation. Therefore, for increasing i7nuci/A*, 
the exponent r will decrease (see Fig.j^b), solid line). 

E. Origins of the discrepancy with Ttheor = 7/6 

The theoretical prediction ED, ^theor = 7/6, is bascd on 
the assumption that the collision rate between two droplets 
of size Si and S 2 can be factorized as described in Eq. ( [^ . 
When one adopts the hypothesis that r takes a universal value, 
this approximation can be conveniently chosen to calculate the 
values of r. When r depends on the microscopic details of the 
dynamics, as observed in Fig.|^ such an assumption must be 
tested for the data at hand. 

To determine the collision rate for our numerical data, we 
count the number of collisions per unit length per unit size 
in two different time intervals Ati = h — to and At 2 = 
h — ti, both in the self-similar regime of the droplets. In 
order to have a better statistics, we average the sets of results 
from five simulations with equivalent initial conditions. These 
numerical data are then compared to the integral of Eq. 
between t and t ^ At. Upon substitution of the expressions 
for n{sR) and N{t) given by Eqs. (|^ andpA]respectively, we 
find 


Anco\\{sl^ s2;t^ At) ^ H {siS2)^^ ^ (t + 1 _ a^^/2 


2-1 


where H 




ymJ \so 


mJ Uso 


Consequently, Eq. (20) implies that 


^<011 = 


Ancoii(sl,s2;i, Ai) 


-ff (SlS2)(^ (sf + 52^"^) (t + At)^^ ^ ^ 


= const. 


(21a) 

(21b) 

( 22 ) 
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FIG. 10. Reduced number density of collisions per unit length An*ou, during two different time intervals Ati = ti—to (left) and At2 = t2 — ti 
(right), with to = 480.5 s, ti = 1035.3 s, t 2 = 1308.9 s (same parameters as in Fig.|^. To a good approximation, the plotted quantities are 
independent of time. However, they show a non-trivial dependence on si and S 2 . 


Hence, the reduced number density of collisions per unit 
length of substrate should be independent of the time t 

and the sizes of the colliding droplets 5i, 52 . In Fig. 10 we plot 
the reduced number density of collisions per unit length of 
substrate An*Qjj during the intervals Afi, At 2 - In particular, 
we use the same set of data of Fig.j^and we take to = 480.5 
s, ti = 1035.3 s , t 2 = 1308.9 s. We find that An*Qjj is indeed 
invariant in time, to a good approximation (see Fig.[TT]). How¬ 
ever, in variance with the assumed expression for the collision 
rate, Eq. ( [^ , it depends on si and 52 . In Fig.[^a) we show a 
cross section of Fig. [T^ along the main diagonal 5i = 52 = 5 , 
both for the first time interval Ati (dashed line) and the sec¬ 
ond time interval At 2 (solid line). In other words, we display 
the reduced number of collisions between droplets of the same 
size. The data present pronounced, very sharp peaks, which 
do not move in time. The first peak, at 5o, corresponds to an 
enhanced probability of having a collision between droplets 
originating at the same moment, upon nucleation in the same 
gap, when they reach the size 5o = (A*/2)^ (see Fig. 1^. The 
second highest peak corresponds to an enhanced probability 
of collision between droplets of size 5 J = A*^. Following 
the mechanism described in Fig.|^ we infer that this peak ac¬ 
counts for collisions between droplets generated at the same 
moment in the same droplet chain, that have already collided 
once. An interpretation along the same lines can be given for 
the peaks appearing at 52 = (2A*)^ and 53 = (4A*)^. With 
our simulation parameters, 5o = 87 /im^ and 5 ^ = 700 //m^, 


52 = 5613 /im^ and 53 = 44900 /im^. The intermediate peaks 
appearing between the main ones can be related to alternative 
merging mechanisms between droplets originated at the same 
moment in the same droplet chain. For example, the peak 
between 5o and 5 ^, corresponds to a size — (3/2A*)^, 
suggesting a mechanism where every second droplet of the 
chain has already merged with its neighbor and the next merg¬ 
ing takes place between droplets of alternated sizes A */2 and 
A*. For colliding droplets larger than 53 , such mechanisms 
are not relevant anymore, thus suggesting the appearance of 
a self-similar area, where the system has lost memory of the 
microscopic details of the nucleation. 

Fig.pH]^b) displays further cross sections of Fig, during 
Ati (dashed lines) and At 2 (solid lines). These cross sections 
are perpendicular to the main diagonal, and they show the 
peaks appearing Fig.[TTJa) at 5o (black lines) and 52 (lit-gray 
lines, green online). The data display a recurring structure, 
with a peak corresponding to the main diagonal, and a dip just 
next to it. Such a structure is present along the whole main 
diagonal, as soon as the droplets are large enough to allow 
for secondary nucleation in the emerging gaps. The peak cor¬ 
responds to the collisions between droplets of the same size. 
The dip surrounding the peak can be explained through the 
concept of self-similarity itself. The global droplet size distri¬ 
bution n(5, t) can be regarded as a superposition of bimodal 
droplet size distributions, with the same shape displayed in 
Fig.[^ but different values of F. If we consider a generic por- 
































































































11 



FIG. 11. (Color online) Cross sections of the data displayed in Fig.[^(a) along the main diagonal and (b) perpendicularly to it, during two 
different time intervals Ati (dashed lines) and At 2 (solid lines), (a) Along the main diagonal, sharply defined maxima appear at sizes s* (with 
z = 0, 1, 2,3). (b) Perpendicularly to the main diagonal, along the lines S 1 S 2 = (so)^ (black lines) and S 1 S 2 = (^ 2 )^ (ht-gray lines, green 

online), there are maxima at the diagonal si/s 2 = 1, dips to its left and right and a plateau for | log^o -si /s 2 1 ^ 3. 


tion of the area occupied by the droplets, where the largest 
size is S*, the rescaled droplet size distribution of such an 
area will present several droplets of similar size (see the 
bump in Fig.[^, from which, the peaks centered on the main 
diagonal in Fig. [T^ originate. Such local droplet size distribu¬ 
tion will also present a gap similar to the one of Fig. from 
which the dips of Fig.[TTJb) originate. 

We conclude that the factorization of the collision rate Acoii 
Eq. ( [20| ), based on the assumption of the distribution of the 
droplets sizes si being uncorrelated to the distribution of the 
sizes of the neighbors S 2 , does not hold for our data. 


V. CONCLUSIONS 

So far, only few numerical 1^ l42l |4^ and experimen¬ 
tal works l(38]| have specifically addressed breath figures on 
a one-dimensional substrate. Performing repeatable and con¬ 
trollable experiments for droplets on a fiber presents technical 
difficulties, such as keeping the temperature of a thin fiber 
constant. Quasi-one-dimensional settings have been realized 
by means of scratches on a plate 13^ , which allowed for a 
better control of the temperature of the substrate. 

From the classical theory of breath figures it is known that 
the size distribution of droplets on a substrate becomes self¬ 
similar after some time that the system evolves; therefore it 
can be described by means of scaling laws, at least in the in¬ 
termediate range of droplets sizes (the poly disperse range). 
In particular, two scaling exponents appear: 0 and r. The 
value of 0 can be inferred from dimensional analysis, while 
the value of the polydispersity exponent r is non-trivial. The 


present work investigated the dependence of this exponent on 
the microscopic details of the system, in order to verify the re¬ 
cent finding ED that r might not be universal, as commonly 
assumed. We developed an event-driven model for three- 
dimensional droplets on a one-dimensional substrate (fiber). 
We included the following details: the growth of the droplets 
by mass deposition from an impinging flow of supersaturated 
vapor, the precursor film between adjacent droplets, the nu- 
cleation process by surface tension driven instability and the 
merging of adjacent droplets. We calculated r, as well as the 
exponents characterizing the decay in time of the porosity and 
the number of droplets. The relations among these three expo¬ 
nents derived in the classical scaling analysis (13116112711311, 
hold also for our model. As expected, we found that the ex¬ 
ponent r does not depend on the impinging mass flow and 
on the monodispersity level of the nucleation process. How¬ 
ever, it does depend on the interaction range 5 of the droplets 
and on the ratio between the nucleation radius -Rnud and the 
spacing A* between the nucleating droplets. In particular, r 
grows with increasing interaction ranges 5 and decreasing ra¬ 
tios i?nuci/A*- Otir results contradicts the expectation that the 
exponent r should be universal. Additionally, the values of 
r that we inferred from our simulations differ by 10 stan¬ 
dard deviations from the theoretical prediction Ttheor = 7/6 
ED Such a prediction was based on the assumption that the 
probabilities of having two neighboring droplets of prescribed 
sizes are uncorrelated. We analyzed the distribution of the 
collisions respect to the sizes of the colliding droplets and we 
showed that this assumption is inaccurate when one keeps into 
account the microscopic details of the system. 

Our observations pose new questions on the non-universal 
nature of the polydispersity exponent r, such as the speciflc 
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mechanisms underlying its parametric dependence as well as 
its quantitative values. These questions have been raised here 
for a setting where the flux impinging on the droplets is pro¬ 
portional to the wetted length on the flber, i.e. a situation 
where the droplet growth is limited by the diffusive mass 
transport to the flber. Alternatively, one could also consider 
a setting where the droplets grow in proportion to their total 
surface area. The classical expectation was that this change 
of the microscopic details of the droplet growth should not 
affect the value of r. In order to test this expectation and to 
address the emerging questions, further numerical simulations 


should be developed, as well as a new theoretical framework 
accounting for the non-universal nature of r. 
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